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ABSTRACT 



We construct the gas-grain chemical network model which includes carbon 
isotopes (^^C and ^^C) with an emphasis on isotopomer-exchange reactions. Tem- 
poral variations of molecular abundances, the carbon isotope ratios (^^CX/^^CX) 
and the isotopomer ratios (^^C^^CX/^^C^'^CX) of CCH and CCS in cold dense 
cloud cores are investigated by numerical calculations. We confirm that the 
isotope ratios of molecules, both in the gas phase and grain surfaces, are signif- 
icantly different depending on whether the molecule is formed from the carbon 
atom (ion) or the CO molecule. Molecules formed from carbon atoms have the 
CX/^^CX ratios greater than the elemental abundance ratio of [^^C/^^C]. On 
the other hand, molecules formed from CO molecules have the CX/^^CX ratios 
smaller than the [^^C/^^C] ratio. We reproduce the observed C^^CH/^^CCH ra- 
tio in TMC-1, if the isotopomer exchange reaction, ^^CCH + H ^ C^^CH -|- H 
-I- 8.1 K, proceeds with the forward rate coefficient kf > 10~^^ cm^ s~^. How- 
ever, the C^^^CS/^^CCS ratio is lower than that observed in TMC-1. We then 
assume the isotopomer exchange reaction catalyzed by the H atom, ^^CCS + 
H ^ C^^CS -|- H -|- 17.4 K. In the model with this reaction, we reproduce the 
observed C^^CS/^^CCS, CCS/C^^^g CCS/^^CCS ratio simultaneously. 

Subject headings: ISM: abundances — ISM: clouds — ISM: molecules 
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1. INTRODUCTION 

The elemental abundance ratio of [^^C/^^C] in the local ISM is 60 (e.g. Langer & 
Penzias (1993); Lucas & Liszt (1998)). The molecular abundance ratio of ^^CX/-'^^CX is 
often assumed to be the same as this elemental ratio. It has been reported, however, that 
the carbon-chain molecules are diluted in ^^C and the ratio varies among molecules; the 
ratio of HCgN/Hi^^CCCN is 79 (Takano et al. 1998), CCS/^^^CCS is 230 ± 130 (Sakai et 
al. 2007) and CCH/i^CCH is larger than 250 (Sakai et al. 2010). These ratios should 
reflect the production pathway of each molecule. In addition, recent observations in TMC-1 
indicated that the abundances of the ^^C isotopomers depend on which carbon atom in a 
molecule is substituted by ^^C. Takano et al. (1998) observed the three ^^C isotopomers 
of HC3N, and found that HCCi^CN is more abundant than HC^^qcn and H^^ccCN, 
which indicates three carbon atoms are not equivalent in HC3N. Takano et al. (1998) 
suggested that the neutral- neutral reaction of C2H2 -|- CN is the main formation path to 
HC3N. Sakai ct al. (2007; 2010) reported the abundance ratios of C^^CS/^^CCS = 4.2 ± 
2.3 and C^^CH/^^CCH = 1.6 ± 0.4. Again, these results indicate two carbon atoms are not 
equivalent in CCS and CCH. Sakai et al. (2007; 2010) suggested that the neutral-neutral 
reactions of CH -|- CS and CH2 + C would significantly contribute to the CCS and CCH 
formation, respectively. These results indicate the neutral-neutral reactions might play an 
important role in production of such carbon-chain molecules. In this way, observations of 
^^C species would be useful to investigate chemistry of organic species. 

Langer et al. (1984) and Langer & Graedel (1989) constructed theoretical models of 
the ^^C chemistry in the gas phase. They showed that carbon bearing species are divided 
into three groups: CO, HCO+ and the "carbon isotope pool". CO is ^^C-rich and "carbon 
isotope pool" is ^^C-poor. HCO"*", on the other hand, can be either ^^C-rich or ^^C-poor. 
Recently Woods & Willacy (2009) revisited carbon isotope fractionation including the 
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adsorption of gaseous molecules onto grain surfaces and grain surface reactions, but they do 
not distinguish ^^C isotopomers. In addition, temporal evolutions in molecular clouds are 
not investigated in detail, because they focused on ^^C species in protoplanetary disks. 

In this paper, we present a detailed chemical model in cold dense cloud cores to 
investigate carbon isotope fractionation and isotopomer fractionation of CCH and CCS. In 
Section 2, we describe our fiducial model. In Section 3, we show the evolutions of molecular 
abundances, isotope ratios and isotopomer ratios in the fiducial model. In Section 4, we 
compare our results with the observations in TMC-1. We also discuss the dependence of 
our results on the assumptions of rate coefficients and initial condition of carbon. Since 
some terms are confusing, we define them here; "carbon isotope ratio" is the abundance 
ratio of the ^^CX/^^CX, and "carbon isotopomer ratio" is the abundance ratio of the 
i2q13q^^i3q12q-^ If a carbon isotope ratio deviates from the elemental [^^C/^^C] ratio, we 
call it " carbon isotope fractionation" . If carbon isotopomer ratio deviates from unity, we 
call it " carbon isotopomer fractionation" . 

2. FIDUCIAL MODEL 

2.1. Carbon isotope fractionation 

There are two processes which fractionate carbon isotopes in molecular clouds: selective 
photo dissociation and isotope-exchange reactions. 

2.1.1. Selective photodissociation 

Carbon monoxide is photodissociated mostly by line absorption. Because of the high 
abundance of CO, self-shielding is efficient especially for ^^CO. Therefore the ^^CO/^^CO 
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ratio is enhanced in the very thin surface region {Ay < 1 niag) (e.g. van Dishoeck et al. 
1988). In this study we consider central regions of molecular cloud cores with Ay ~ 10 
mag, and hence we do not consider this process in the following sections. 

2.1.2. Fractionation through isotope-exchange reactions 

There are small zero-point vibrational energy differences among the isotopically distinct 
species. In a cloud of low temperature {T ~ 10 K) the isotope fractionation occurs through 
isotope-exchange reactions. Watson et al. (1976) pointed out that the carbon isotope 
fractionation in interstellar molecules mainly occurs as a result of the reaction: 

where AE, the zero-point vibrational energy difference between the reactants and products, 
is 35 K. The rate coefficient of this reaction was measured by Watson et al. (1976) at 300 
K and by Smith & Adams (1980) at 80 - 510 K. Langer et al. (1984) extrapolated these 
experimental data below 80 K assuming the following relation between the forward reaction 

rate coefficient kf and the backward reaction rate coefficient k^. 

kr = kf exp{—AE/kBT), 

where /cb is Boltzmann coefficient. At 10 K, the forward rate coefficient is 1.34 x 10~^ cm^ 

s~^ and the backward rate coefficient is 4.0 x 10^^^ cm'^ s~^. Smith & Adams (1980) also 
measured the rate coefficient of the following proton transfer reaction at 80 - 510 K: 

H^^co^ + ^^CO ^ ^^CO + H^^co^ + AE, (2) 

where the vibrational zero-point energy difference is AE = 9 K. Langer ct al. (1984) 
estimated the reaction rate coefficient at 10 K: kf = 6.5 x lO"^'^ cm^ s~^ and kr — 2.7 x 
10-^° cm=^ s-^ 
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2.2. Carbon isotopomer fractionation 

Sakai et al. (2007; 2010) discussed the cause of carbon isotopomer fractionation of CCS 
and CCH based on their observational results in TMC-1. They pointed out two possibihties: 
(i) the formation pathways of the species and (ii) the exchange of the ^^C position after 
formation of molecules by isotopomer-exchange reactions. Specifically they proposed the 
reactions in Table 1. 



2.2.1. Fractionation through formation paths 

If the observed fractionation is due to the formation paths, the following neutral-neutral 
reaction (e.g. Turner et al. 2000) could significantly contribute to the CCH formation: 

C + CH2 ^ CCH + H. 

Two carbon atoms are not equivalent in this neutral-neutral reaction: 

^^C + ^^CH2 ^ ^^C^^CH + H, (3) 

i^C + ^^CHs ^ ^3Qi2Qfj ^ fj_ 

If the isotope ratios of the carbon atom and CH2 are different, the isotopomer ratio of CCH 
deviates from unity. 

Similarly, the observed isotopomer fractionation of CCS could be due to the 
neutral-neutral reaction: 

CH + CS ^ CCS + H, 
in which two carbon atoms are not equivalent: 

^^CH + ^^CS ^ ^^C^^CS + H, (5) 

^^CH + ^^CS ^ ^^C^^CS + H. (6) 
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In addition, the isotopomer fractionation of CCH might propagate to CCS through the 
following reactions (Yamada et al. 2002): 

CCH + S ^ CCS + H, (7) 

and 

S+ + CCH ^ CCS+ + H, (8a) 

CCS+ + H2 ^ HCCS+ + H, (8b) 

HCCS+ + e ^ CCS + H. (8c) 

In the reaction (7), the sulfur atom will attack the end carbon atom of CCH, which has 
the unpaired electron, to form the CCS structure (Yamada et al. 2002). In the case of the 
sulfur ion, the reaction may occur in a similar way. So C^^CS is made from ^^CCH, and 
^^CCS is made from C^^CH through these processes. If these processes are efficient, the 
C^'^CS/^'^CCS ratio should be less than unity, given the C^'^CH/^'^CCH ratio is greater than 
unity. This is inconsistent with the observation: both the isotopomer ratios of CCH and 
CCS are greater than unity. 



2.2.2. Isotopomer- exchange reactions 

Fractionation by isotopomer-exchange reactions is also possible: 

^=^C^^CH + H ^ ^^C^^CH + H + A£;, (9) 

where A£J, the zero-point vibrational energy difference of ^^CCH and C^^CH, is 8.1 K 
(R. Tarroni, private communication; Tarroni & Carter 2003). As a candidate of the 

isotopomer-exchange reaction of CCS, on the other hand, Sakai et al. (2010) considered the 
following reaction: 

^=^c^^cs + s ^ ^^c^^cs + s + A£;, (10) 
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where A£', the zero-point vibrational energy difference of ^^CCS and C^^CS, is 17.4 K. 
According to the energy surface calculations, these exchange reactions, (9) and (10), do not 
have the activation barrier (Y. Osamura private communication). At low temperature (T ~ 
10 K), the isotopomer ratios of CCH and CCS can be greater than unity. 

2.3. Reaction network 

We construct a chemical reaction network with carbon isotopes and isotopomers. 
The network is based on the gas-grain model developed by the Ohio State University 
astro chemistry group (Garrod & Herbst 2006), and is expanded to include isotopes and 
isotopomers by the three procedures described below. While Garrod & Herbst model 
consists of 655 species and 6309 reactions, our reaction network consists of 1203 gas and 
grain species, and 21116 gas phase and grain surface reactions. 

Firstly, we include ^^C species into the reaction network statistically. For reactions 
involving ^^C species, we assume that the total rate coefficient is the same as that of the 
corresponding reaction of the main isotope. We distinguish isotopomers only in species with 
less than three carbons to save computational time. For example, ^^C^^C^^CH, ^^C^^C^^CH 
and are different species, but i^Ci^Ci^Ci^CH, i2ci3ci2ci2CH, i^ci^ci^ciscH 

and ^^C^^C^^C^^CH are treated as the same species, ^^Cg^CH. We checked molecular 
structures, so that symmetric molecules, such as C2H2, do not have isotopomers. We do not 
consider multiple ^^C species for simplicity. The statistical expansion increases the number 
of the species and reactions by about two and three fold, respectively. For example, the 
reaction 

CCH + O ^ CH + CO, 
with rate coefficients of k is replaced by the following five reactions with rates k, k/2, k/2, 
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A;/2, and /c/2, respectively: 

i^C^^CH + O -)■ ^^CH + ^^CO, 
i^C^^CH + O ^ ^^CH + ^^CO, 
^^C^^CH + O ^ ^^CH + i^CO, 
i^C^^cH + O ^ ^^CH + i3co. 

In this way, statistical expansion is the main cause of increasing the number of species and 
reactions. "Statistically" means no fractionation occurs. In other words, all the isotope 
ratios are equal to the elemental [^^C/^^C] ratio and all the isotopomer ratios are unity. 

Secondly, we include the gas phase reactions in Table 1 referring to the discussion in 
Section 2.2. The rate coefficients of these reactions are not measured in the laboratory. So 
we assume that the rate coefficients are 1 x 10~^° cm^ s~^ for neutral-neutral reactions and 
1 X 10^^ cm^ s~^ for ion-neutral reactions. The dependence of our results on these rate 
coefficients is discussed in section 4.3. 

Finally, we modify the branching ratios of proton transfer reactions and electron 
recombination reactions which break the C-H bond, such as 

H^^C^^CS^ + e^^^C^^CS + H, 

fji2ci3cs+ + e^ 12^13^8 + H, 

In the statistically expanded network, the rate coefficients of these reactions are set to k', 
k'/2, and k' /2, respectively, where k' is the rate coefficient of the recombination of the ^^C 
species. Then these reactions tend to make the isotopomer ratio of CCS unity, even if 
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the isotopomer ratio of HCCS''' are fractionated. But it would be more likely that these 
reactions do not rearrange the structure of C-C-S. So we set the rate coefficient to be k' , k' , 
and zero, respectively. 

Table 2 lists the elemental abundances in our model; it is known as "low metal" values 
because of their strong depletions of species heavier than oxygen (e.g. Graedel et al. 1982). 
Initially the species are assumed to be atoms or atomic ions except for hydrogen, which is 
in molecular form. We assume the elemental isotope ratio, [^^C/^^C], of 60. If it is higher 
or lower than 60, the resultant molecular isotope ratios are scaled by [^^C/^^C]/60. We 
confirmed in the range of values, [^^C/^^C] = 30 - 77, the deviation from this scaling is less 
than 4 %. 

We use the cosmic ray ionization rate of 1.3 x 10"^'^ s~^, which is commonly used for 
chemical model simulations for dense molecular clouds (e.g. Terzieva & Herbst 1998). We 
assume that the sticking probability of colliding gaseous species onto grains is unity, and 
the adsorbed species can diffuse on the grain surfaces by thermal hopping. We adopt the 
same adsorption energies and non-thermal desorption rates as Garrod & Herbst (2006). As 
a physical model, we assume a static cloud core. In our models, the physical parameters are 
set to T = 10 K, riHa = 5 x 10"^ and Ay = 10 mag, referring to the observation of TMC-1 
(Pratap et al. 1997), unless otherwise stated. In section 3, we describe the results of our 
fiducial model. 
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3. RESULT 
3.1. Isotope ratios 

3.1.1. Isotope ratios of gas species 

In this paper, we define 'abundance' as the fractional abundance with respect to H 
nuclei. Figure la shows the evolutions of major carbon-bearing species. Initially the carbon 
exists as carbon ion, which recombines to form carbon atoms in 10^ yr. After a few 10^ yr 
major form of carbon is mainly in carbon monoxide, which depletes onto the grain surfaces 
after 5 x 10^ yr. 

Figure lb shows the evolutions of the isotope ratios of the carbon-bearing species. Our 
results agree with Langer et al. (1984); carbon bearing species are divided into three groups: 
CO, HCO"'" and the "carbon isotope pool". Although we include gas-grain interactions, it 
does not affect the isotope ratios in the gas phase before the depletion of CO {t < 10^ yr). 
Adsorption rates are inversely proportional to square root of the mass of species, so that 
the rate of ^^C species and ^^C species are different. But this difference is small (< 4 %) 
and do not affect significantly. It should be noted, however, that when the CO depletion 
becomes significant, the reaction (1) becomes less efficient and isotope ratios of "carbon 
isotope pool" species are lowered compared to the model without gas-grain interactions. 
Since the thermal desorption is ineffective at 10 K, grain surface reactions do not affect 
isotope ratios of gaseous species. Then the isotope fractionation of the gaseous species is 
mostly determined by the gas phase reactions and the depletion of CO. Before a few 10^ yr 
the abundance ratio of CO/^^CO is significantly smaller than 60 because of the reaction (1). 
As the abundance of CO increases with time and CO becomes the dominant carbon-bearing 
species, the isotope ratio gets close to the elemental abundance ratio of [^^C/^^C] = 60. For 
example, the CO/^^CO ratio is 54 at 10^ yr. In contrast, the C'^/^^C+ ratio is larger than 
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60 because of the reaction (1). Many carbon species are produced by reactions starting 
from C+ or C, which is produced by recombination of C+. Therefore the isotope ratios of 
these "carbon isotope pool" species are similar to the C"'"/^'^C"'" ratio before 10^ yr and the 
C/^^C ratio after 10^ yr, although actual values vary among species depending on their 
production pathways. For example, the CN/^^CN ratio is 83, while CS/^^CS ratio is 109 at 
10"^ yr. The evolution of the HCO"''/H^^CO+ ratio is similar to those of "carbon isotope 
pool" at t < 10^ yr, but the ratio is lowered by the reaction (2) after a few 10^ yr; the 
HC0+/H13C0+ ratio is 45 at 10^ yr. 

To investigate the density dependence of the isotope ratios, we also performed 
calculations with = 5 x 10^ and 5 x 10^ cm~^. The evolution is in general faster at 
higher densities. Since we are interested in carbon-bearing species, we compare isotope 
ratios of these molecules at the time of peak carbon-chain abundance; 1 x 10^, 2 x 10^ and 
4 X 10^ yr for densities of 5 x 10^, 5 x 10^ and 5 x 10^ cm~^, respectively. We found that 
the isotope ratios do not significantly depend on density. For example, the CN/^^CN ratio 
at the selected time is 83, 81 and 80, and the CS/^^CS ratio is 83, 92 and 108 for densities 
of 5 X 10^, 5 x 10"^ and 5 x 10^ cm~^, respectively. The isotope ratios of other species in 
each density model are listed in Table 3. 



3.1.2. Isotope ratios of ice mantle species 

Figure 2 (upper panels) shows the evolution of the major ice mantle species. H2CO and 
CH3OH are the most abundant organic species in ice mantles. They are formed from the 
carbon atom at ^ < 2 x 10^ yr and from CO at later time. 

Figure 2 (lower panels) shows the evolution of the isotope ratios of the ice mantle 
species. Carbon atoms and carbon monoxides in the gas phase are adsorbed and form more 
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complex species by grain surface reactions. Before 10^ yr, the isotope ratio of carbon atoms 
on grain surfaces is almost the same as in the gas phase. Carbon monoxide, on the other 
hand, is directly absorbed on grain surfaces and is also formed by the grain surface reactions 
of the carbon-bearing species. Then the isotope ratio of CO on grain surfaces is shghtly 
different from that in the gas phase. The isotope ratios of ice mantle species depend mainly 
on whether the species are formed from the carbon atom or the CO molecule; the isotope 
ratio is larger than 60 if the species are formed from C atom, while the ratio is smaller 
than 60 if the species are formed from CO. For example, CH3OH is formed from the carbon 
atom and its isotope ratio is 81 at t = 2 x 10^ yr. After 2 x 10^ yr, CH3OH are formed 
from the CO molecule and the isotope ratio decreases to 53 at 10® yr, which is almost the 
same as the isotope ratio of CO. HCOOH is an exception. It is formed from HCO"*" in the 
gas phase and adsorbed onto the grain surfaces. The isotope ratio of HCOOH is similar to 
that of HCO^ in the gas phase. Isotope ratios of ice species in higher (5 x 10^ cm^'^) and 
lower (5 x 10^ cm~^) density models are listed in Table 4. 

Boogert et al. (2000; 2002) found that the solid C02/^"^C02 ratio is similar to the solid 
CO/^^CO ratio towards NGC 7538 IRS 9, and suggested that CO2 is formed from CO on 
grain surfaces, rather than from the carbon atom in the gas phase. In our model, CO2 is 
formed in the gas- phase via the reaction of an O atom with HCO, which is formed from 
carbon atom, and adsorbed onto grain surfaces. Then the isotope ratio of CO2 is similar 
to that of the carbon atom. However if we use the modified rates method (e.g. Ruffle and 
Herbst 2000), which artificially slows down the surface diffusion rate of the hydrogen atom, 
CO2 is formed from CO on grain surfaces via the reaction of HCO -t- O (Aikawa et al. 
2005). In this case, the isotope ratio of CO2 ice at t > 2 xlO^ yr should be closer to that of 
CO ice just like H2CO and CH3OH. 

It should be noted that the observed isotope ratio of CO2 ice varies with lines of sight. 
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The solid CO/^^CO and COs/^^COs ratios towards NGC 7538 IRS 9 are actually 71 ± 15 
and 80 ± 11, respectively (Boogert et al. 2000; 2002). The isotope ratio of CO2 is similar to 
or slightly higher than that of CO. On the other hand, towards NGC 2264 IRSl, Langer & 
Penzias (1990) found the gas C^^O/^^^iso ratio is 56 ± 3, and Gibb et al (2004) found the 
sohd 002/^^002 ratio is 131 ± 21. The isotope ratio of CO2 ice is higher than that of CO 
gas by a factor of 2.3, assuming the carbon isotope fractionation of ^^O species is similar to 
that of ^^O species (e.g. Langer et al. 1984). In our model, in which CO2 is mainly formed 
in the gas phase and adsorbed onto grain surfaces, such high isotope ratios of CO2 ice 
relative to that of CO ice (gas) is reahzed at 10^ yr in the model with rin^ = 5 x 10^ cm~^. 
These variations may suggest that interstellar CO2 ice is formed from both in the gas phase 
and on grain surfaces, and contributions of the two paths varies among clouds. 

3.2. Isotopomer ratios 

3.2.1. CCH 

Figure 3a shows the temporal evolution of the assorted species relevant to the 
isotopomer fractionation of CCH and CCS. Figure 3b shows the C^^CH/^^CCH ratio. 
The isotopomer ratio is significantly enhanced in the fiducial model (solid line). In order 
to investigate whether the isotopomer fractionation of CCH is caused by (i) formation 
reactions or (ii) isotopomer-exchange reactions, we performed a calculation without the 
exchange reaction (9) (dashed line in Figure 3b). It is clear that isotopomer fractionation is 
mostly due to the exchange reaction. 

Figure 4 shows the major formation and destruction paths of CCH isotopomers. The 
reactions (3) and (4) are the dominant formation reactions until lO'' yr. The isotope ratio 
of CH2 is larger than that of carbon atoms, and hence the isotopomer ratio is lowered. 
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After 10^ yr, on the other hand, CCH is mainly produced from the destruction of larger 
species (e.g. CsH^) by the oxygen atom and electron recombinations of symmetrical species 
(e.g. C2H3 ). These reactions tend to make the isotopomer ratio of CCH unity, since they 
form the same amount of the CCH isotopomer: C^^CH and ^^CCH. Therefore, without the 
exchange reaction (9), the isotopomer ratio becomes almost unity in the model. 

In our model, we do not distinguish isotopomers in species with more than four 
carbons. We checked if this simphfication affects our results. We performed calculations 

without these large carbon chains (more than 4 carbon atom) and the exchange reaction 
(9), and confirmed the isotopomer ratio of CCH is as low as the dashed line in Figure 3b. 

In the model with the exchange reaction, the isotopomer ratio increases with time in 
a few 10^ yr as the electron abundance decreases. Then, the ratio gradually decreases as 
the H atom abundance decreases by adsorption and subsequent grain surface reactions. At 
several 10'' yr, the ratio starts to increase again by the adsorption of oxygen atoms. After 3 
X 10^ yr, it reaches nearly the equilibrium ratio of 2.2 {kf/kr). 

Since the H atom abundance is an important factor for the efficiency of the exchange 
reaction (9), we briefly discuss the formation and destruction of H atoms in dense cloud 
cores (Figure 3a). Before 10^ yr, the primary sources of the H atom are the reactions 
of hydrocarbons, such as the recombination of CH^ and the reaction of CH + O, and 
the destruction is by adsorption on grain surfaces followed by subsequent hydrogenation 
reactions, such as H + OH. At t > 10'' yr, H atoms are mainly formed by cosmic ray 
induced photodissociation of hydrogen molecules. At this stage the absolute abundance of 
the H atom is ~1 cm""^ independent of the gas density. We also checked the dependence of 
our results on the initial H atom abundance; if the initial abundance is lowered by an order 
of magnitude, the isotopomer ratio varies only slightly. 

Figure 5 shows density dependence of the C^^CH/'^^CCH ratio. As the density 
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increases, both the abundances of electron and H atom decrease. While the electron 
abundance is proportional to n^^ ^, the H atom abundance is proportional to n^^. Then the 
peak isotopomer ratio decreases with density. 

The isotopomer fractionation also affects the isotope ratios, i.e. the CCH/^^CCH and 
CCH/C^^CH ratios (Figure 6). The reaction (9) enhances ^^C in C^^CH, while it dilutes 
i3q I^CCH. For comparison, the dashed line in Figure 6 shows the isotope ratios of CCH 
in the model without the isotopomer fractionation: the model without the reaction (3), (4), 
(5), (6), (9) and (10). 



3.2.2. CCS 

Figure 7 shows the temporal evolution of the C^^CS/^^CCS ratio. The isotopomer 
ratio of CCS is lower than unity in the model without the exchange reaction (10) (dashed 
line). The reactions (5) and (6) are not the main formation reactions of C^^CS and ^^CCS 
in our model. Instead, CCS is mainly formed by the electron recombination of larger ions, 
such as HC2S"'" (Figure 4). These reactions tend to make the isotopomer ratio of CCS unity. 
The reaction (7) and (8) decreases the isotopomer ratio of CCS, because the isotopomer 
ratio of CCH is larger than unity. In the fiducial model (solid line) , the isotopomer ratio is 
lower than unity except for very limited periods. The isotopomer ratio reaches a maximum 
of about 1.2 by the reaction (10) at 3 x 10^ yr. Figure 8 shows the density dependence 
of the C^^CS/^^CCS ratio. Although the peak isotopomer ratio increases with density, 
the maximum ratio is 1.8, which is much lower than observed. Unlike CCH, the exchange 
reaction (10) is inefficient, because the sulfur atom is less abundant than the hydrogen atom 
by three orders of magnitude. 
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4. DISCUSSION 

4.1. Compcirisons with observations of TMC-1 

The gas density in TMC-1 is 10^ - 10^ cm'^^ (Pratap et al. 1997). In the model with 
riHa = 5 X 10^ cm~^, the abundances of carbon-chain molecules, such as HC3N, reach their 
maxima at 2 x 10^ yr. Since carbon-chain molecules are abundant in TMC-1, we compare 
the obtained molecular abundaces and isotope and isotopomer ratios in the fiducial model 
at 2 X 10^ yr to the observed values in TMC-1. 

Pratap et al. (1997) made maps along the TMC-1 ridge, and determine the column 
densities of C^^O and CCH. The molecular abundance of CCH in TMC-1 is 5.7 x 10"^° - 
4.0 x 10~^, assuming a C^^0/H2 ratio of 1.7 x 10"'' (Frerking et al. 1982). In our model, 
we obtain the abundance of CCH to be 3.0 x 10~^, which is in reasonable agreement with 
the observed value. 

In the fiducial model, we obtain the CCH/^^qch, CCH/C^^CH and C^^CH/i^CCH 
ratios to be 99, 63 and 1.6, respectively, at 2 x 10^ yr. The observed CCH/C^^CH and 
CCH/i^CCH ratios are higher than 170 and 250, respectively (Sakai et al. 2010). The 
average isotope ratio, which is calculated by 

CCH 

^ ^ C13CH + 13CCH' 

is higher than 200, and is three times higher than our model. Additional isotope-exchange 
reactions, which are currently unknown, may be needed to reproduce the observed values. 
On the other hand, the observed C^^CH/^^CCH ratio is 1.6 ± 0.4 (Sakai et al. 2010), which 
is in reasonable agreement with our model. Therefore the exchange reaction (9) can account 
for the isotopomer fractionation of CCH. 

We obtain the CCS/^^CCS, CCS/C^^qs and C^^cs/iSQCS ratios to be 84, 76 and 1.1, 
respectively, in the fiducial model at 2 x 10^ yr. The observed CCS/^^CCS, CCS/C^^^CS 
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and C^^CS/^^CCS ratios are 230 ± 130, 54 ± 5, and 4.2 ± 2.3, respectively (Sakai et al. 
2007). The average isotope ratio is 87, which is in good agreement with our model results. 
But our model underestimates the isotopomer fractionation. We discuss this problem in the 
next section. 



4.2. Exchange reactions of CCS 

The isotopomer-exchange reaction (10) is not efficient enough to account for the 
observed isotopomer ratio of CCS, because of the low abundance of sulfur atoms. Since 
the reactant is not necessarily the sulfur atom for the exchange reaction, we consider a 
possibility of the following neutral-neutral reaction: 

^^C^^CS + H ^ ^=^C^^CS + H + AS, (11) 

where AE' is 17.4 K. This can be regarded as a catalytic reaction by the hydrogen atom. 
According to the quantum chemical calculation by Yamada et al. (2002), the reaction 
(11) does not have the activation barrier. It could be more effective than the reaction 
(10), because the hydrogen atom is more abundant than the sulfur atom. Since the rate 
coefficients of the reaction (11) has not been measured or calculated, we simply assume that 
the forward rate coefficient is 1 x 10"^" cm^ s~^, which is a typical coUisional rate coefficient 
between neutral species. Figure 9a shows the temporal evolution of the isotopomer ratio of 
CCS in the model with the exchange reaction (11). The ratio is significantly enhanced in 
comparison with the model in Figures 7 and 8, and reaches nearly the equilibrium ratio 
of 5.7 (kf/kr) at a few 10^ yr. After several 10^ yr, CCS abundance decreases by about 
two orders of magnitude (Figure 3a). The abundance of S-bearing molecular ions, such as 
HCCS+, decrease only by an order of magnitude, because abundance is enhanced by the 
CO depletion. Then the exchange reaction (11), whose efficiency is proportional to the CCS 
abundance, becomes less effective than ion-molecule reactions, and the isotopomer ratio 
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decreases to 4.1. Our model with the reaction (11) can reasonably explain the observed 
isotopomer ratio. 

The isotopomer fractionation also affects the isotope ratios, CCS/^^CCS and 
CCS/C^^^CS (Figure 9b). is significantly diluted in ^^CCS, and the CCS/^^^CCS ratio is 
higher than 200 in several lO"' yr. On the other hand, the CCS/C^^CS ratio is lower than 
60. At 2 X 10^ yr, the CCS/^^qcs and CCS/C^^CS ratios are 263 and 47, respectively, 
which are in reasonable agreement with the observed isotope ratios. Hence the reaction 
(11) is a possible culprit for the isotopomer fractionation of CCS. 

4.3. Dependences on rate coefficients of isotopomer-exchange reactions 

In the previous sections we assumed the forward reaction rate coefficients of the 
isotopomer-exchange reactions are all 1 x 10~^° cm^ s~^, which is a typical coUisional 
rate coefficient between neutral species. However it is not obvious whether the actual 
exchange reactions proceed at the collisional rate. For example, the rate coefficient of the 
similar type of reactions, such as HNC3 + H ^ HC3N + H, is assumed to be 1 x 10~^^ 
cm~^ s~^ in our model (e.g. Garrod & Herbst 2006). In particular, reaction (11) requires 
the rearrangement of the order of the carbon nuclei in the chain through the transient 
cyclic isomers (Yamada et al. 2002). Although the reaction is exothermic, the efficiency of 
reaction (11) might be low. Hence, we investigate the dependence of the isotopomer ratios 
on the rate coefficients of reactions (9) and (11). Figure 10 shows the temporal evolution 
of the isotopomer ratios of (a) CCH and (b) CCS. The forward reaction rate coefficients 
are 1 x 10~^° cm^ s~^ (sohd hne), 1 x 10~^^ cm^ s~^ (dashed hne) and 1 x 10~^^ cm^ s~^ 
(dotted line). We can see the forward reaction rate coefficients should be larger than 1 x 
10~^^ cm^ s~"^ to account for the observed isotopomer ratios by the exchange reactions (9) 
and (11). Laboratory measurements, as well as more detailed theoretical calculations of 
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these isotopomer-exchange reactions, are highly desirable. 

4.4. Initial condition of the Ccirbon 

So far, we have assumed that all the carbon initially exists as carbon ions. This 
assumption is the most favorable for the carbon isotope fractionation, because the 
isotope-exchange reaction (1) becomes efficient. Here, we investigate the dependencies of 
isotope ratios on the initial condition of the carbon. 

We performed calculations with two different initial conditions; (1) a half of the carbon 
exists as atomic carbon, and (2) all the carbon is in atomic form. In the former case, the 
evolution of isotope ratios in the gas phase is very similar to the fiducial model, because 
the abundances of major carbon-bearing species are almost the same as the fiducial model 
after 10^ yr. Figure 11 shows the evolutions of (a) molecular abundances and (b) isotope 
ratios in the gas phase for the latter case: all the carbon is initially atomic carbon. The 
molecular abundances and isotope ratios are different from the fiducial model before 10^ 
yr, but they approach to the fiducial model after a few 10^ yr. The abundances of C, C"*" 
and CO are similar to those in the fiducial model after 10^ yr. The isotope ratio of CO is 
nearly 60 before 10^ yr, since the abundance of C"*" is small. After 10^ yr the ratio shghtly 
decreases: ^^CO/^^CO = 56 at 10^ yr. On the other hand, the isotope ratio of C"*" is as high 
(> 60) as in the fiducial model, since the abundance of CO is smaller than in the fiducial 
model only by an order of magnitude. The isotope ratio of the carbon atom is nearly 60 
before 10^ yr. After 10^ yr the ratio becomes greater than 60 because of the neutralization 
of the carbon ion. "Carbon isotope pool" species are mostly produced from the carbon 
atoms, then their isotope ratios are similar to that of the carbon atom. In this way, the 
carbon isotope fractionation occurs by the isotope-exchange reactions after 10^ yr, even if 
all the carbon is initially in the carbon atom. These results suggest isotope fractionation 
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in early times does not much affect tlie isotope ratios after 10^ yr. In order to account for 
the isotope fractionation at the period of peak carbon chain abundances ~ 2 x 10^ yr), 
however, some fraction of carbon should be carbon ion initially. 

4.5. Summciry 

We constructed a gas-grain chemical network model that includes carbon isotopes 
and isotopomers in order to investigate the evolution of molecular abundances, the carbon 
isotope ratios and the isotopomer ratios of CCH and CCS in cold dense cloud cores. The 
principal results are follows. 

1. We confirm that the isotope ratios of molecules, both in the gas phase and on grain 
surfaces, mostly depend on whether the species is formed from the carbon atom (ion) or the 
CO molecule; the isotope ratio is larger than the elemental abundance ratio of [^^C/^^C] 

if the species is formed from the carbon atom, while the ratio is smaller if the species is 
formed from the CO molecule. 

2. There are two possible processes of isotopomer fractionation: (i) the formation 
pathways of the species and (ii) the exchange of the ^"^C position by isotopomer-exchange 
reactions. We reproduce the observed isotopomer ratio C^^CH/^^CCH by considering the 
isotopomer-exchange reaction: ^^CCH -|- H — >■ C^^CH -|- H -|- 8.1 K, if this reaction proceeds 
with the forward rate coefficient kf > 10~^^ cm^ s~^. The isotope ratio is underestimated 
in our model, which is left for future studies. 

3. We reproduce the observed C^^^CS/^^^CCS, CCS/C^^^CS and CCS/^^^CCS ratios 
simultaneously by considering the isotopomer-exchange reaction: ^^CCS -|- H — ;> C^'^CS + 
H -I- 17.4 K, if this reaction proceeds with the forward rate coefficient kf > 10~^^ cm^ s~^. 

4. In conclusion, isotopomer fractionation of CCH and CCS can be due to isotopomer- 
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exchange reactions. Laboratory measurements and detailed quantum calculations of these 
isotopomer exchange reactions are highly desired. 
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10^ 10^ 10"^ 10^ 10^ 

Time [yr] 



Fig. 1. — Temporal variation of (a) molecular abundances and (b) isotope ratios in the gas 
phase in the fiducial model. 
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10^ 10^ 10* 10^ 10** 10^ 10^ 10* 10^ 10^ 

Time [yr] Time [yr] 



Fig. 2. — Temporal variation of molecular abundances (upper panels) and isotope ratios 
(lower panels) of ice mantle species in the fiducial model. 
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Fig. 3. — (a) Temporal variation of the abundances of CCH, CCS, H atom, O atom and 
electron, (b) Temporal variation of the C^^CH/^'^CCH ratio in the fiducial model (solid line) 
and in the model without the exchange reaction (8) (dashed line). 
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Fig. 4. — Production and destruction pathways of CCH and CCS. 
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Fig. 5. — Isotopomer ratio of CCH in the model with = 5x10^ (dashed hne), 5x10"^ 
(sohd hne), and 5x10^ cm~^ (dotted hne). 
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10^ 10^ 10"* 10^ 10® 

Time [yr] 



Fig. 6.— Temporal variation of isotope ratios of CCH: the CCH/i^QCfj g^^d CCH/C^^^CH 

ratios. The sohd hncs show ratios in the fiducial model, whereas the dashed line shows the 
ratio in the model without isotopomer fractionation. 
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Fig. 8. — Isotopomer ratio of CCS in the model with 71^.2 — 5x10^ (dashed hne), 5x10^ 
(sohd hne), and 5x10^ cm~^ (dotted hne). 
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io^ 10^ 10"^ 10^ 10^ 

Time [yr] 



Fig. 9. — Temporal variation of (a) isotopomer ratio, C^^CS/^^CCS and (b) isotope ratios 
of CCS. The solid lines show the ratios in the fiducial model, while the dashed line shows 
the isotope ratio of CCS in the model without the isotopomer fractionation. 
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10^ 10^ 10"^ 10^ 10^ 



Time [yr] 

Fig. 10.— Temporal variation of the (a) C^^cH/i^CCH ratio, and (b) Ci^CS/i^CCS ratio. 
Forward rate coefficients of the reaction (9) and (13) are 1 x 10~^° cm^s~^ (sohd fine), 1 x 

10"^^ (dashed Unc), and 1 x 10"^^ cm^s^^ (dotted hne). 
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Fig. 11. — Same as Figure 1, but all the carbon is in carbon atoms initially. 
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Table 1: Rate coefficients and zero point energy differences of tlie reactions added to the 
statistically expanded network 



Reaction 


k {cm^ s ^) ^ 


AE (K) 


CH + CS ^ CCS + H 


i.oo(-io)i 




CH + ^ C^^CS + H 


I.OO(-IO) 




i^CH + CS ^ i^ccs + H 


I.OO(-IO) 




CCS+ + H2 ^ HCCS+ + H 


1.00(-9) 




C13CS+ + H2 ^ HC13CS+ + H 


1.00(-9) 




13CCS+ + H2 ^ H13CCS+ + H 


1.00(-9) 




13C+ + CO ^ i^co + c+ 


1.34(-9) 


35 


HCO+ + i^CO ^ CO + H13C0+ 


6.50(-10) 


9 


13CCH + H ^ C^^CH + H 


I.OO(-IO) 


8.1 


i^CCS + S c^^cs + s 


I.OO(-IO) 


17.4 



Note. — (a) For exchange reactions k are the rate coefficients of forward reactions, (b) x{—y) represents 
X X lO-^'. 
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Table 2: Initial abundances with respect to hydrogen nuclei 



species Abundance 

H2 5.00(-l) 

H 5.00(-5) 

He 9.75(-l) 

1^0+ 7.86(-5) 

1.31(-6) 

N 2.47(-5) 

1.80(-4) 

Si+ 2.74(-9) 

S+ 9.14(-8) 

Fe+ 2.74(-9) 

Na+ 2.25(-9) 

Mg+ 1.09(-8) 

C1+ 1.00(-9) 

P+ 2.16(-10) 
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Table 3: Isotope ratios and fractional abundances of gaseous species at 10*^, 10^ and 10^ yr 
in the model with riRa = 5 x 10^, 5 x 10^ and 5 x 10^ cm~^ 



Species 


riH, (cm ^) 


10" yr 
rii/nu 




10^ yr 
rii/nu 


/ 


10^ yr 
ni/riH 


12CX/13CX 


CO 


5x10^ 


6.0(-6) 


14 


2.8(-5) 


43 


2.5(-5) 


57 




5x10^ 


1.7(-5) 


34 


1.6(-5) 


54 


1.3(-6) 


50 




5x10^ 


1.2(-5) 


51 


2.0(-7) 


48 


1.2(-7) 


49 


c 


5x10^ 


6.5(-5) 


79 


2.6(-5) 


80 


8.1 (-9) 


154 




5x10^ 


4.0(-5) 


79 


1.6(-8) 


84 


3.0(-8) 


77 




5x10^ 


7.0(-7) 


74 


9.7(-10) 


65 


5.2(-10) 


66 


HCO+ 


5x10^ 


3.4(-ll) 


88 


2.3(-9) 


51 


5.7(-9) 


52 




5x10'' 


4.9(-10) 


51 


1.9(-9) 


45 


5.0(-10) 


50 




SxlO'"" 


2.7(-10) 


42 


9.8(-ll) 


48 


7.3(-ll) 


49 


H2CO 


5x10^ 


5.9(-10) 


158 


1.7(-9) 


74 


3.3(-9) 


76 




5x10'' 


2.9(-10) 


76 


3.2(-10) 


74 


2.4(-9) 


68 




5x10^ 


5.6(-ll) 


73 


1.6(-10) 


56 


1.3(-10) 


56 


cs 


5x10^ 




150 


2.6(-9) 


83 


1.9(-10) 


109 




5x10^ 


1.4(-8) 


109 


3.0(-10) 


77 


5.1(-10) 


75 




5x10^ 


1.9(-9) 


87 


4.4(-ll) 


72 


4.8(-ll) 


73 


CN 


5x10^ 


7.0(-9) 


130 


8.0(-9) 


83 


3.3(-9) 


182 




5x10^ 


4.6(-9) 


83 


7.6(-10) 


84 


7.7(-9) 


77 




5x10^ 


9.8(-10) 


77 


6.9(-10) 


67 


3.7(-10) 


69 


CH 


5x10^ 


9.0(-9) 


127 


8.0(-9) 


78 


8.7(-10) 


97 




5x10^ 


5.0(-9) 


78 


5.0(-10) 


77 


1.5(-9) 


75 




5x10^ 


2.6(-10) 


74 


l.l(-lO) 


71 


1.3(-10) 


72 




X iU 




1 no 




78 
1 


O.Z[-LL ) 


00 




5x10" 


3.5(-9) 


79 


1.2(-10) 


75 


2.5(-10) 


76 




5x10-'^ 


6.7(-10) 


73 


9.1(-11) 


72 


l.l(-lO) 


73 


HCN 


5x10^ 


3.6(-10) 


151 


1.4(-8) 


93 


5.3(-9) 


180 




5x10" 


4.6(-8) 


85 


1.5(-9) 


96 


l.l(-8) 


77 




5x10^ 


4.1(-8) 


77 


9.6(-10) 


68 


4.1(-10) 


69 


HNC 


5x10^ 


3.0(-10) 


172 


1.3(-8) 


96 


5.3(-9) 


186 




5x10" 


3.6(-8) 


88 


1.3(-9) 


101 


l.l(-8) 


77 




5x10^ 


4.0(-8) 


77 


9.2(-10) 


67 


3.8(-10) 


69 
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Table 3: continue 







10^ yr 




10^ yr 




10^ yr 




Species 


riHs (cm-^) 


rii/nu 


i2CX/"CX 


rii/riH 






i^CX/i^CX 


HC3N 


5x10^ 


4.8(-13) 


123 


9.4(-10) 


85 


1.2(-11) 


99 




5x10^ 


1.8(-9) 


76 


1.4(-10) 


73 


2.8(-ll) 


75 




5x10^ 


9.5(-10) 


71 


2.9(-13) 


73 


9.4(-14) 


71 


HC5N 


5x10^ 


7.2(-14) 


149 


7.2(-ll) 


83 


1.1(-12) 


88 




5x10^ 


4.3(-10) 


76 


9.6(-12) 


71 


2.2(-12) 


74 




5x10^ 


9.6(-ll) 


68 


1.6(-14) 


73 


4.3(-15) 


89 


CCH 


5x10^ 


6.8(-9) 


136 


8.6(-9) 


80 


7.2(-lG) 


107 




5x10^ 


3.3(-9) 


78 


2.2(-9) 


76 


8.3(-10) 


76 




5x10^ 


8.2(-10) 


73 


2.7(-ll) 


71 


2.4(-ll) 


70 


CCS ^ 


5x10^ 


2.0(-10) 


140 


1.0(-9) 


81 


3.9(-10) 


97 




5x10* 


1.3(-9) 


83 


2.4(-10) 


74 


2.3(-ll) 


75 




5x10^ 


6.2(-10) 


80 


2.9(-13) 


81 


3.4(-13) 


79 


C3S 


5x10^ 


7.0(-12) 


149 


I.O(-IO) 


82 


3.4(-12) 


86 




5x10^ 


l.l(-lO) 


82 


3.2(-ll) 


71 


2.5(-12) 


74 




5x10^ 


1.2(-10) 


73 


1.8(-14) 


79 


1.6(-14) 


81 



Note. — 



(a) Isotope ratios for CCH and CCS are average isotope ratios of isotopomers. 
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Table 4: Isotope ratios and fractional abundances of ice mantle species at 10^, 10^ and 10^ 
yr in the model with riHa = 5 x 10^, 5 x 10^ and 5 x 10^ cm~^. 







iU yr 




1 fib ,rT* 

iU yr 




1 n6 try 

iU yr 




Species 


riHs (cm ) 


ni/riii 


ISnv /13r<Y 
OA./ l^A 


ni/nn 


OA./ OA 


ni/uH 


OA/ OA 


(^u ice 


X iU 


2.5(-7) 


00 


Z.o[-D j 


/ 


2.1(-5j 


00 




OX iU 


2.7(-bj 


00 


O.Z(-D ) 


oU 




oi 




OX iU 


z. 0(^-0 j 


40 


4.1)^-0 j 




7/" 


A 


t^xiaUri ice 


X iU 


Z.i>[-0 ) 


(0 


7.2(-8j 


oU 


l.b(-b) 


KA 
04 




X iU 




Q1 

oi 


6.ti(-() 


oi 


z.z^-D ) 


04 




5x10^ 


7.4(-7) 


82 


6.2(-7) 


81 


1.3(-6) 


52 


HCOOH ice 


5x10^ 


1.1(-18) 


129 


3.4(-12) 


53 


2.5(-10) 


59 




5x10* 


5.3(-13) 


53 


2.9(-10) 


46 


3.0(-10) 


49 




5x10^ 


4.8(-ll) 


45 


1.5(-10) 


44 


l.l(-lO) 


44 


CO2 ice 


5x10^ 


7.7(-10) 


74 


6.1(-9) 


81 


9.0(-8) 


78 




5x10* 


9.3(-9) 


81 


2.2(-7) 


78 


1.7(-7) 


76 




5x10^ 


3.2(-7) 


79 


5.2(-7) 


78 


2.8(-7) 


77 


H2CO ice 


5x10^ 


6.5(-8) 


75 


2.4(-7) 


60 


1.4(-5) 


54 




5x10* 


4.2(-7) 


81 


9.1(-7) 


80 


2.0(-5) 


53 




5x10^ 


1.3(-6) 


80 


2.6(-6) 


61 


1.5(-5) 


50 


CH4 ice 


5x10^ 


7.2(-7) 


78 


8.2(-6) 


82 


1.2(-5) 


79 




5x10^ 


8.1(-6) 


81 


1.8(-5) 


80 


2.0(-5) 


78 




5x10^ 


2.2(-5) 


81 


2.2(-5) 


81 


2.3(-5) 


80 



